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ABSTRACT 

We use the self-consistent model technique developed by Martm-Manjon et al. 
(2008) that combines the chemical evolution with stellar population synthesis and 
photo-ionization codes, to study the star formation scenarios capable of reproducing 
the observed properties of star-forming galaxies. 

The comparison of our model results with a database of Hn galaxies shows that 
the observed spectra and colors of the present burst and the older underlying popu- 
lation are reproduced by models in a bursting scenario with star formation efficiency 
involving close to 20 per cent of the total mass of gas, and inter-burst times longer than 
100 Myr, and more probably around 1 Gyr. Other modes like gasping and continuous 
star formation are not favored. 

Key words: galaxies: evolution - galaxies: star formation -galaxies: starburst - 
galaxies: ISM - ISM: HII regions 
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1 INTRODUCTION 

Hn galaxies can be considered as the strong emission line 
subset of the Blue Compact Dwarf galaxies (BCD). Their 
optical emission is dominated by strong and narrow emis- 
sion lines produced by the interstellar gas ionized by a 
young and luminous cluster. Their integrated optical spec- 
tra are indistinguishable from a normal Hn region. Hn 
galaxies are gas rich and metal poor. These facts had 
led to the early proposal that these systems may be very 
young per haps undergoing their first burst of star forma- 
tion (SF) (|Sargent fc Searlej . ll970l ). However, many authors 
have found clear evidences of the existence of a low surface 
brightness, old, non-ionizing stellar po pulation in the major- 
ity of t hese galaxies, as explained by iMartm-Manion et al.l 
l|2008al .here in after MMDT) and references therein. The 
existence of an old population plus the relative paucity of 
Hn galaxies with extreme low metal content are consistent 
with a scenario in which Hn galaxies are suffering at present 
a starburst, instead of being old systems with a slow evolu- 
tion. Therefore we use the term starburst to denote a violent 
episode of star formation where a large number of massive 
stars have been formed in a small volume of space and over a 
time scale of a few million years. Galaxies called Starbursts 
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Galaxies are a very heterogeneous category, including from 
BCDs to ULIRGS. We will handle here only Hn galaxies, 
those dwarf galaxies with active star formation whose ef- 
fects actually dominates their UV-optical emission. In these 
galaxies the starburst is considered a phase, a process of 
strong star formation in terms of intensity and duration. 
Nevertheless, the fundamental details of how the star for- 
mation history proceeds in Hn galaxies is still an unresolved 
problem. 

Until now, three basic evolutionary scenarios have been 
postulated: 

• Bursting star formation: the stars form in short 
but intense episodes separated by long quiescent peri- 
ods of very low or null activity (|Davies fc Phillippsl . 1 19881 : 
iBradamante et~aL 1. 1 19981 ) 

• Gasping star formation: The star formation takes 
places as long episodes of SF of m oderate in t ensity - 
separated by short quiescent periods dTosi et all Il99ll : 
lAparicio fc Gallartl , Il995l : iRecchi fc Henslerl |2004 ) ■ 

• Continuous star formation: The process of stellar 
formation is continuous and of low intensity d uring the 
galax y life, with superimposed sporadic bursts (|Legrandl . 

ibood) . 

The low metallicities, lack of dust and the optical colors 
argue in favour of a bursting star forma tion with long qui- 
escent periods (e.g. iMarconi et all 1 1994 ). However, inactive 
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star forming periods seem to rarely occur, and inter-burst 
stages characterized by low-levels of star formati on are per- 
haps more plausible than its complete cessation (|Lee et al.l . 
|2004 ). If this is the case, a continuous low level state of star 
formation would dominate during the vast majority of the 
time during which many or perhaps most of the stars in 
BCD galaxies would be formed. 

Chemical evolution models with gasping or burst- 
ing star formation modes seem t o be appropriate to de 



scribe local dw a rf irre g ulars (dlrr) dRecchi fc Hensler , 2004 
Marconi et all, 1 19951 ; iTosi et all Il99ll ; iGallagher et al 



199^), but the chemical evidence alone cannot differenti- 



ate between gasping and bursting scenarios. The gasping 
scenario can be viewed as a transition betwe en the continu- 
ous s tar formation and the bursting scenario l|Marconi et al.l . 
1 1994 ). The most remarkable difference between gasping and 
bursting scenarios is that short episodes of star formation 
enrich the ISM in a time scale of few tens of Myr while 
long-lasting episodes enrich gradually the ISM in a longer 
time scale, and any further episode of SF does not leave an 
appreciable imprint on the chemical evolution. The abun- 
dances, which increase rapidly in bursting models, change 
more slowly and during a longer time in gaspi ng models, but 
with older ages in their stellar populations |Recchi et all 
|2001| . |2002| ). This picture, however, does not entirely exclude 
the possibility that star formation histories appear to be 
composed by burst cycles when examined with a high time 
resolution, but looks constant when averaged over longer 
time scales. It does not exclude either the possibility that the 
SF propagates through the galaxy, taking place in indepen- 
dent luminous short lived Hn regions: "the short timescales 
associated to starbursts in dwarf galaxies may be understood 
as Bickering ev ents, small component s of a larger starburst 
in the galaxy" (|McQuinn et all 120091 ) 

These observations of stellar populations suggest that 
it is most likely that the star formation in dwarf galax- 
ies is sporadic, separated by millions to billions of years, 
even in isolated systems, as most of Hn are. If the SF of 
Hn galaxies is dominated by intermittent starburst episodes 
separated by long inactive time spans, quiescent blue dwarf 
galaxies without the dominant starburst and showing similar 
properties to those with such intense star formation events, 
should be relatively common. Support to this framework 
comes from the fact that blue compact dwarfs (BCD), with 
emission lines on average much weaker than these of Hn 
gala xies, are more common than t hose ones. Following this 
idea lSanchez Almeida et al.l (|2008l ) obtain a relationship be- 
tween the duration time of the starburst phase and the qui- 
escent periods, based on the number of objects of each type 
of their sample. They find that, if the duration of the burst 
phase is 10 Myr, the time in quiescence must be at least 
0.27 Gyr, implying several mayor starburst episodes along 
the life of a BCD galaxy. Therefore, the diversity in proper- 
ties exhibited by dwarfs may be at least partially the con- 
sequence of observing them at different times in their star 
formation cycles, and the frequency distributions of galaxies 
in the various phases could correspond to the time spent in 
each of them. 

It would be interesting to see which is the most prob- 
able star formation scenario given by the ab initio cos- 
mological hydrodynamical simulations. However this is dif- 
ficult since the dwarf galaxies produced in the existing 



works are usually poor gas, dwarf elliptical (dE) or dwarf 
spheroidals, (dSp) objects JValcke. de Riicke. fc Deionghd . 
120081 ; iRevaz et all [20091; ISawala et alll201Ch. Other s imula- 
tions, such as iPelupessv. van der Werf. fc Ickd \2004 ) used 
fixed initial conditions for stellar and gas masses, and this 
way the m ass to halo ma s s rati o is not a result of the sim- 
ulation. In IStinson et al.l (|2009l ) simulations, the Ultravio- 
let (UV) background is not included, which probably con- 
tributes to a hi gh star formation efficiency and to lar g e fina l 
stellar masses. iMashchenko. Wadslev. fc Couchmanl (|2008l ) 
follow the evolution of a galaxy with a halo mass of 10 9 M 
but only until z=5. iGovernato et al] (|2010l ) have calculated 
some complete hydrodynamical simulations, in which bary- 
onic processes, as gas cooling, heating from the cosmic UV 
field, star formation and supernova-driven gas heating, are 
included with sufficient spatial resolution, with clumps as 
low as 10 M© resolved. The created dwarf galaxies with- 
out bulges are similar to the observed dwarf irregular galax- 
ies. However, besides the authors do not show the star for- 
mation histories, these galaxies show a maximum star for- 
mation rate (SFR) of O.25M yr _1 and a present SFR of 
0.01 MQyr^ 1 , much lower than observational values for Hn 
galaxies, and not comp arable to these ob jects that we will 
try to model here. Only iNagamind (|2010l ) has recently pre- 
sented some star formation histories, obtained from cos- 
mological hydrodynamical simulations for dwarf galaxies, 
with enough spatial resolution for gravitational masses of 
10 9 Mq. Their star formation histories seem to be sporadic, 
and the stars continue forming sporadically even at late 
times. However, the simulations are still having some prob- 
lems to reproduce the adequate number of galaxies and 
the author c l aims that more work is still necessary. As 
ISawala et al.l (|201ll ) explains, all dwarf galaxies ( 10 10 Mq) 
formed in the current hydrodynamical simulations are more 
than an order of magnitude more luminous than expected 
for these masses. In any case they are much more massive 
than BCDs and Hn galaxies that we want modeling. 

For what refers to other models existing in the liter- 
ature, there are not models for the study of BCD or HII 
galaxies in the way we want to use them. A number of 
models have computed purely chemi cal evolution models 
for BCDs or dwarf i rregu l ar galaxies |Chiosi fc Matteuccil . 
19821 ; iMarconi et al.l.ll994l;lRecchi et alll2002l;l Recchi et all 



200a IShiet all 120061 ; iRecchi fc Henslerl , 120071 , and many 
others). Most of them assume that the star formation occurs 
in bursts and include the effects of galactic winds and/or gas 
infall. However, they limit the study to the evolution of N 
and O abundances and/or t o the lumi nosity-metallicity rela - 
tion. iMouhcine fc Continil (|2002l ) and I Vazquez et all (|2003l ) 
used the inf ormation coming fr om the chemical evolution 
models from ICarigi et al.l (|2002T l to perform the next step 
and combine chemical and spectral evolution for irregular 
galaxies. These models, however, exclude the early stages of 
evolution, i.e. during the nebular phase when most massive 
stars do minate the energy output budget. | Kruger et al.l 
Jl99ll); iLindner. Fritze-v. Alvensleben. fc Frickd (1 19991) ; 



IVazdekis et all (|l997l ) or Ivan den Hoek et al.1 (|2000l ) com- 
pute chemical and photometric evolution models in a way 
more or less consistent for starburst, spiral, early types 
or low surface brightness galaxies, respectively, any of 
them applied to Hn galaxies. Other works, in turn, are 
focused on the ionized gas properties and make models 
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using Single Stellar Populations (SSPs) spectral energy 
distributions (SEDs) of a given metallicity for applying 
a photo-ionization code to obtain the emission lines and 
studying diagnostic diagrams an d /or abundances ob- 
taine d by empirical calibrations dStasiriska, fc Leitherer 
19961: IStasinska. Schacrcr, & Leitherer, l200ll : 



20061 ; 



Stasinska fc Izotovl 120031; IMov et all , l200ll : iDopita et al 



IMartm-Manion et all 



20 id ). 



Most of them ignore 



the star formation history of the galaxy and not take into 
account the underlying stellar population, so these studies 
are only valid for the study of the current stellar generations 
of the galaxy. Summarizing a code combining chemical, 
evolutionary synthesis and photo-ionization models has not 
been applied for the analysis of BCDs and Hn galaxies, with 
the exception of the ours (MMDT), but only to e lliptical 
massive galaxies (jBressan. Chiosi. fc Fagotto! . 1994 ). 

The aim of this work is to critically analyze the possibil- 
ity of dwarf galaxies undergo recurrent phases of star forma- 
tion by comparing the predictions of our self-consistent evo- 
lutionary models with the most evolution sensitive observed 
parameters, i.e. stellar continuum colors, emission line equiv- 
alent widths and chemical composition in a large sample of 
BCD galaxies. 

We have used an updated version of our models from 
MMDT, developed to predict the main characteristics of Hn 
galaxies, such as the emission lines, continuum color, con- 
tinuum plus the contribution of emission lines colors, equiv- 
alent widths and chemical abundances. Our code uses the 
chemical evolution model results for the computation of the 
SEDs, which are, in turn, used as the ionizing sources for 
a photo-ionization code. In MMDT we studied the viability 
of this tool, showing its adequacy in the reproduction of the 
data defined by the stellar populations (colors, equivalent 
width of absorption lines or spectral indices, spectral energy 
distributions..), which define the time evolution, as well as 
the characteristics of the gas phase (emission lines, equiv- 
alent widths, elemental abundances, gas densities), which 
defines the present time state of the galaxy. This is done in 
a self-consistent way, that is, using the same assumptions 
regarding stellar evolution, model stellar atmospheres and 
nucleosynthesis, and using a realistic age-metallicity rela- 
tion. 

One of the most important results obtained from our 
previous work was that observational data are reproduced 
only if the mass involved in the last burst is much smaller 
than the mass of the old underlying stellar population. Now, 
once checked the possibilities of our tool, we try to wide the 
number of models, changing the star formation scenarios to 
see if data are reproduced with more or less success for some 
of them. 

It is a widely assumed idea that starburst galaxies are 
strongly affected by gas infall and outflow and that chemi- 
cal abundances and star formation cannot modeled without 
these issues. However, the existence of outflows in dwarf star- 
burst galaxies is by no means a settled issue; while winds 
able to escape the galaxy have been found in two proto- 
typical starburst galaxies like NGC1569 and NGC1705, re- 
cent works a bo ut ma ss loss, for example Ivan Evmeren et al.l 
(|2007l . l2009l lbl l2010l ), find that ionized and neutral hydro- 
gen expansion velocities measured are, in all cases, too low 
to allow the gas to escape from the gravitational potential 
of their studied galaxies some of which, like NGC2363 or 



NGC4861, c an also be considered as prototypical. In fact, 
according to lBomans et al.l (|2007l ):" While the observational 
support for the presence of galactic winds in massive galax- 
ies and gas-rich mergers is quite strong, the case for galactic 
winds in dwarf galaxies is much weaker" . 

The need of winds is also related with the low oxygen 
abundance found for dwarf galaxies for their gas fractions, 
and, in this case, both infall and outflow can help. However, 
galaxies loosing a large fraction of their gas should show red 
colors, contrary to what is observ ed. In fact, a good fra ction 
of the galaxies in the sample of Ivan Zee et alj (|2006l ) can 
be reproduced with closed box models. Furthermore, it is 
necessary to take into account the N/O vs O /H relationship 
for dw a rf galaxies. Selective wind s models as lMarconi et al] 
l|l994 ) ; iBradamante et all (|l99ST ) have recognized problems 
in reproducing the observed N/O ratio s, and the work by 
lLarsen. Sommer- Larsen, fc Pagel (|200ll) demonstrates that 
the observed N/O ratios, as well as their dispersion, can be 
reproduced with closed box models while s elective winds can 
be ruled out. Actually. iMolla et alj (|2006l ) showed that it is 
mainly the combination of the different time scales in nitro- 
gen production by stars of different masses and the galaxy 
star formation histories as proceeding from different star 
formation efficiencies, what establishes the N/O ratio. 

Therefore, since there is enough proof to conclude that 
infall and/or outflows are not necessary to reproduce the 
general trends of Hn galaxies, galactic supernova-driven 
winds are not included in our models. This scenario is also 
in agr eement with models from lTassis. Kravtsov. fc Gnedir] 
(]2008i ) , who find that scaling relations of dwarf galaxies may 
be reproduced by simulated galaxies without supernova- 
driven outflows. 

For the present work, new and updated theo- 
retical codes have been applied for the computa- 
tion of the models. The main changes correspond 

to the evolutionary synthe sis models^ now from 

IMolla. Garcia- Vargas fc Bressanl (l2009l.her e inafte r 



MGVB09) instead iGarcfa-Vargas.Bressan. fc Dfazl (|l995l) . 
As explained before, to guarantee the self-consistency of the 
approach we have used the same assumptions in the stellar 
evolution, model stellar atmospheres and nucleosynthesis 
parts of the code. The use of these new models allows to 
be more consistent than in our previous work MMDT, 
since there the IMF used in the chemical evolution models 
dFerrini et all. ll990T ) and in evolutionary synthesis models 



(|Salpeterl . ll955l ) was not the same. The new models MGVB 
allow to us use the same IMF for chemical and evolutionary 
synthesis (photometrical) calculations. Moreover, with 
the p revious models from IGarcfa-Vargas.Bressan. fc Dfazl 
(|l995l ) we had also a problem for the lowest metallicities 
stellar populations since models for the youngest ages of low 
metallicities were not available. Now we have models for 
the same number of ages for all metallicities. To include the 
lowest metallicities of the youngest ages change appreciably 
the colors of the stellar populations and allows to calculate 
emission l ines for very low m e tallicit y regions, as it is 
shown in IMartm-Manion et alj (|2010l ). The contribution 
of the emission lines is quite different than expected from 
the simple extrapolation from other metallicities, also 
modifying the final results of colors. 

The data shown in the graphs have been extracted 
mainly from two main sources. First, the compilation from 



© 2011 RAS, MNRAS 000, [TV?? 



4 Martin-Manjon et al. 



iHovos fc Di'azl (|2006h which provides emission line measure- 
ments, corrected for extinction, published for local Hn galax- 
ies. The sample comprises 450 objects and constitutes a 
large sample of local Hn galaxies with good-quality spec- 
troscopic data. The sample is rather inhomogeneous in na- 
ture, since the data proceed from different instrumental se- 
tups, observing conditions and reduction procedures, but 
have been analysed in a uniform way. Data for these sam- 
ple objects include the emission line intensities of: [Oil] A A 
3727,29 A [OIII]AA 4959,5007 A and [NII]AA 6548,84 A all 
of them relative to H/3, and the equivalent widths of the 
[Oil] and [OIII] emission lines - EW([OII]) EW([OIII]) - 
, and the H/3 line, EW(H/3). As a second source, we have 
used the metal poor galaxy data from t he Data Release 3 o f 
Sloan Digital Sky Survey, taken from llzotov et al.l (|2006h . 
The Sloan Digital Sky Survey (|York et all |2000| ) consti- 
tutes a large data base of galaxies with well defined se- 
lection crite ria and observed in a homogeneous way. The 
SDSS DR3 (|Abazaiian et all . l2005h provides spectra in the 
wavelength range fro m 3800 to 930 A fo r ~ 530000 galax- 
ies, quasars and stars. ITzotov et all (|2006h extracted ~ 2700 
spectra of non-active galaxies with the [OIII]A4363 A emis- 
sion detected above la level. This initial sample was fur- 
ther restricted to the objects with an observed flux in the 
H/3 emission line larger than 10~ 14 erg s _1 cm" 2 and for 
which accurate abundances could be derived. They have 
also excluded all galaxies with both [OIII]A4959/H/3 < 0.7 
and [OII]A3727/H,g > 1.0. Applying all these selection cri- 
teria, they obtain a sample of ~ 310 SDSS objects. Data for 
these sample objects include the emission line intensities of: 
[OIII]A4959,5007 A and [NIIJA6584 relative to H/3 and the 
equivalent width of H/3. They also include the intensity of 
the [Oil] AA 3727,29 A emission line for the lowest redshift 
objects. 

In the next section an explanation of the theoretical 
star-bursting models is made. The section 3 presents the 
results, the meaning of each input parameter and the im- 
plication of their variation over the results of the models. 
We will discuss the influence of these parameters in the re- 
production of the observable characteristics of star-forming 
galaxies, their impact over possible star-formation scenarios 
and the connection among different type of dwarf galaxies 
in section 4. Finally, a summary and the conclusions of this 
work are presented in section 5. All tables of the models will 
be available in electronic format. 



2 THEORETICAL MODELS 

The viability of our model technique to reproduce the ob- 
servable characteristics of Hn galaxies was discussed in 
MMDT: In a first stage the chemical evolution model is 
computed by using certain input parameters to obtain the 
star formation history and the gas and stars chemical abun- 
dances. Secondly, the evolutionary population synthesis code 
is applied to obtain the SEDs corresponding to each time 
step of the chemical evolution. Finally, the ionizing part of 
the SED and the resulting abundances are used as an input 
to the photo-ionization code to compute the time evolution 
of the emission lines of the ionized gas. 



2.1 Chemical evolution 

The chemical evolution is computed with a simplified ver- 
sion of the classica l mult ip hase chemical evolutio n mode l 
from iFerrini et all (119941 ); iMolla. Ferrini. fc Diaz! (|l99fj ): 
iMoIla fc Di'azl (|2005l 'l. In our version of the code there is 
only one region (without the two zones halo and disk) 
with a given mass of gas which form stars (in this case we 
do not consider the phase of formation of molecular clouds 
as in the classical multiphase models). We have run models 
considering the star formation as a set of successive bursts 
followed by quiescent periods in a region with a total mass 
of gas of 10 8 Mq . In each burst a given amount of gas is 
consumed to form stars, and this process comes defined by 
the star formation efficiency. The code solves the chemical 
evolution equations to obtain, in each time step, the abun- 
dances of 15 elements: H, D, 3 He, 4 He, C, 13 C, O, N, Ne, 
Mg, Si, S, Ca, Fe, and nr (where nr are the isotopes of the 
neutron rich elements, synthesized from 12 C, 13 C, 14 N and 
16 O inside the CO core). The stellar yields are those from 
IWooslev fc Weaver! (1 19951) for massive sta rs, M > 8 Mq , and 
those from lGavilan. Buell. fc Mollal (|2005h for low and inter- 
medi ate mass stars. The supernova la yields used proceed 
from llwamoto et al. Jl99£ ). The initial mass function (IMF) 
is taken from lFerrini et al (|l99Cf ) with a r ange between .15 
and 100 Mq. M ore recent IMFs such as iKroupal (|2002T ) or 
IChabrierl (|2003l ) there exist but, for chemical evolution mod- 
els, it is necessary to use a combination IMF+stellar yields 
sets calibrated with Milky Way Galaxy data, that is, able to 
reproduce the elemental abundances and gas and star den- 
sities as observed. This calibration is not still done for those 
IMFs, although an update of models taking these IMFs and 
some new stellar yield sets will be shown in a next future 
(Molla et al. in preparation). We have taken time steps of 
<5t~O.7Myi0 from the initial time, t = 0, up to the final one, 
t = 13.2 Gyr. At each time step, the star formation rate and 
the mass in each phase - low mass, massive stars and rem- 
nants, total mass in stars created, and mass of gas- are also 
obtained. 



2.2 Evolutionary synthesis 

The SEDS are taken from the PopStar evolutionary syn- 
thesis mod els by MGVB09. Isochrones are an update from 
those from iBressan.Granato. fc Silval (|l998! ) for 6 different 
metallicities: Z = 0.0001, 0.0004, 0.004, 0.008, 0.02 and 
0.05. The very low metallicity model of Z = 0.0001 had 
not been included before in similar works. The age coverage 
is from logt =5.00 to 10.30 with a variable time resolution 
of A(logi) = 0.01 in the y oungest stellar age s. We have 
used the Ferrini IMF results (I Ferrini et all Il990t ) with mass 
limits between 0.15 and 100 Mq in order to avoid any incon- 
sistency with the chemical evolution code, which also uses 
this IMF. 

To calculate the SED integrated over the whole history 
of the galaxy, SEDs of SSPs with the corresponding metal- 



1 It implies that there is no infall of gas over a disc, as in the 
spiral and irregular galaxies models 

2 The time step is chosen to include the fastest evolutionary 
phases of the most massive stars 
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licity and age must be convolved with the star formation 
history (SFH): 

Lx(t) = f Sx(r,Z(t'))<f(t')dt' (1) 
Jo 

where r = t — t' is the age of the stellar population created in 
a time t and S\ is the SED for each SSP of age r and metal- 
licity Z reached in that time t'. A SED from the SSP library, 
S\ , must be assigned to each time step according to its corre- 
sponding age and metallicity. Taking into account the SFH, 
ty(t), and the age-metallicity relation, Z(t), obtained from 
the chemical evolution model, we know the age and metallic- 
ity assigned to each time step or stellar generation. For each 
one of them we have chosen the SED closest in age among 
those available in the grid of PopStar. However, in our mod- 
els the metallicity changes continuously while the available 
SEDs of the library are computed only for 6 possible values. 
Therefore, we have interpolated logarithmically between the 
two SSP of the same age r and closest in metallicities to Z(t) 
to obtain the corresponding S\(t, Z(t')). The final result is 
the total luminosity at each wavelength A, La (erg.s— l.A -1 ), 
corresponding to the whole stellar population, including the 
ionizing continuum proceeding from the last formed stellar 
population. 



2.3 Photo-ionization calculations 

To compute the photo- ionization models, we have used 
CLOUDY (version c06.02 iFerland et ail 119981 ^1. The gas is 
assumed to be spherically symmetric around a point source 
of radiation and the pressure or density in the gas is im- 
posed by external conditions. This allows a plane-parallel 
geometry treatment in which the gas may be regarded as a 
thin shell. A closed geometry has been taken for the calcu- 
lations. All the photons which escape from the illuminated 
face of the cloud towards the star, go on to strike the other 
side of the nebula, ensuring the case B of recombination and 
the approximation on the spot. The number of ionizing pho- 
tons, Q(H), striking the illuminated face of the cloud, have 
been calculated directly from the total resulting SED of the 
models. 

We derive the radius of the modelled region from the 
mechanical energy from massive stars with strong winds 
(taken from MGVB09) , instead of using the radius required 
to maintain a constant density of stell ar mass through the 
successive bursts, as we did in MMDT. ICastor et all (jl975l ) 
demonstrated that an early-type star with a strong stellar 
wind can blow out a large cavity or bubble in the surround- 
ing gas, if it is assumed to be compressed into a thin spher- 
ical shell. The wind-driven shell begins to evolve with an 
initial phase of free expansion followed by an adiabatic ex- 
pansion phase, and then the material collapses into a thin, 
cold shell as a result of radiative cooling. At this stage the 
gas traps the ionization front and the radiative phase begins. 



3 In order to maintain the consistency with MMDT, we have 
used the same Cloudy version for the present work. The differ- 
ences between the used version and the laters are new features 
and new possibilities for the modelization, however, there are not 
significative differences in our model results. 



In this phase the ionizing photons are absorbed and the re- 
gion cools via emission in the Balmer lines. In this process, 
the radius of the outer shock, R s , evolves as: 

R s = 1.6(Z, mec /n) 1/5 i 3/5 pc (2) 

where L mec is the total injected mechanical energy (SN and 
stellar winds) per unit time in units of 10 36 ergs s _1 , n is the 
interstellar medium density in units of cm' 3 , and t is the 
age of the shell in units of 10 4 yr. Since we use the value of 
the energy at each time-step, this radius represents the in- 
stantaneous size of the region obtained by adding the winds 
and SNe from the previous age to this age. Therefore L mec 
is this energy divided by the time step. Then, the ionized 
gas is assumed to be located in a thin spherical shell at that 
distance R s from the ionizing source. This approach has the 
advantage of eliminating the ionization parameter as a free 
variable in the models since now it is computed from the 
physical parameters of the evolving young cluster. 

The ionized gas abundances are assumed to be those 
reached at the end of the starburst previous to the cur- 
rent one. Fifteen element abundances have been introduced 
in the photo-ionization code: He, C, N, O, Ne, Na, Mg, 
Al, Si, S, Ar, Ca, Fe and Ni, obtained from the chemical 
evolution model, except for Na, Ar and Ni which are not 
computed in the mo del and are scaled to the solar ratio 
jAsplund et al.l . l2005h . The models assume that the nebula 
is ionization bounded and no dust has been included in the 
chemical evolution models neither in the photo-ionization 
calculations. However dust grains mixed with the ionized 
gas have been partially taken into account, since we have 
included the depleti on in refractor y elem ents (Si, Fe, CA, 
Si, Mg) taken from iGarnett et al.l (1 19951 ). The grains can 
affect the absorption of the UV photons and decrease the 
electronic temperature. The density has been assumed con- 
stant for simplicity and equal to 1 00 cm" 3 , whi c h is a ppro- 
priate for modeling Hn galaxie s (lHagele et al. I. l2008h and 
large circumnucl ear Hn regions (|Garcfa- Vargas et all Il997l : 
iDi'az et alJ . l2007h . frequently found around the nuclei of star- 
bursts and AGNs. Although the constant density hypothesis 
is probably not realistic, it can be considered representative 
when the integrated spectrum of the nebula is analyzed. 

The shape of the ionizing continuum is defined by the 
pair of values (i/(Ryd), \oguL v ) (Eq. [1} obtained from the 
ionizing spectrum given directly by the evolutionary synthe- 
sis code. 

2.4 Input parameters 

Taking into account the simplified version used for our chem- 
ical evolution models, the free input parameters related with 
the infall time-scale and the molecular cloud formation dis- 
appear. The total mass of the initial gas is within the region 
from the initial time. The star formation rate is zero except 
during the bursts. We need to define the intensity of these 
bursts, and the time elapsed between them. 

Each model is therefore characterized by three input 
parameters: 

• 1 - The Initial efficiency (e): It is the amount of gas 
consumed to form stars in the first burst of star formation, 
that is, *(t) = = eM g 

We present here models computed with 2 values of e: 
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(i) High efficiency models: The first burst of star forma- 
tion involves 33 per cent of the total initial mass of gas 
(33xl0 6 M Q ), 

(ii) Low efficiency models: The first burst of star forma- 
tion involves 10 per cent of the total initial mass of gas 
(lOxlO 6 M ). 

• 2 - Attenuation: The star formation efficiency of the 
following bursts is attenuated in two different ways: 

(i) By a factor which changes with the number of the 
burst, n, according to the expression: 

*n = (-) ■ *0 

n 

This attenuation type corresponds to a soft attenuation 
model (hereinafter SAM). 

(ii) By a constant factor, fc'™ -1 ', according to the expres- 
sion: 

corresponding to a high attenuation model (hereinafter 
HAM). In this last expression k is the attenuation factor 
which, in this caseQ, takes values from to 1. The lower 
the k value, the stronger the attenuation, the burst being 
less efficient each time; the higher the k value, the weaker 
the attenuation, and then the stronger the bursts. We only 
show models with k = 0.65 

• 3 - Time between bursts (At): Every burst takes 
place instantaneously and it is followed by quiet periods, 
whose duration can change. For this work we have taken 
At= 1.3 Gyr for the inter-burst time, that is, one burst every 
1.3 Gyr as the generic case, although we will also show some 
models with At = 0.1 Gyr and At — 0.05 Gyr in order to 
compare the different results. 

The models we show here are a selection of the ones cal- 
culated and described more widely in Martin-Manjon, Molla 
& Lopez-Sanchez, 2012 (in preparation). In that work, the 
results of 20 models are compared with data for generic and 
particular low mass dwarf galaxies. We have selected the 
most relevant of them in order to show the effects of the 
input parameter variation on the resulting evolutionary his- 
tory of each galaxy. The input parameters of these six mod- 
els are shown in Table [TJ 



Table 1. Input parameters for the theoretical models. In column 
1 it is denned the type of model according to the attenuation: 
soft (SAM) or high attenuation (HAM); The second column gives 
the parameter k which defines the attenuation; column 3 is the 
star formation initial efficiency e, and column 4 shows the time 
between burst At in Gyr. 



Num. Attenuation k e At 

of model type Gyr 

1 SAM - 010 L30 

2 SAM 0.33 1.30 

3 HAM 0.65 0.10 1.30 

4 HAM 0.65 0.33 1.30 

5 HAM 0.65 0.33 0.10 

6 HAM 0.65 0.33 0.05 
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Figure 1. SFR of the models of the table [JJ as labelled. In a) we 
show the SAM (1 and 2), in b) two HAMs (3 and 4), with the 
2 same efficiencies as in a). For these 4 models the time between 
bursts is At = 1.3 Gyr. In c) we show two models (5 and 6) 
with shorter time between bursts, At = 0.1 and 0.05 Gyr. The 
dashed (green) lines define the upper and l ower limits to the S FR 
estimated for BCD and/or Hn galaxies bv lHovos et al. 



3 RESULTS 

We will show here the main results obtained with the se- 
lected models. The efficiency determines the initial star 
formation rate and the initial metallicity of the gas. We have 
plotted in Fig. Q]the SFR for the 6 models. The first burst 
is strong in all panels, while the subsequent ones are less 
intense due to the decrease of the available gas to form stars 
and the attenuation. In a) and b) the last time computed 
is 13 Gyr while in c) and d) the final time is around 1.3 
and 0.6 Gyr respectively. In all panels we show the observa- 
tional limits as dashed (green) lines. Taking these limits into 

4 It may be higher than 1 in models with increasing efficiency; 
we have computed some of these models, however, their results 
are in clear disagreement with most observations. 



account, the dwarf galaxies may suffer 11 bursts for SAM, 
or 8-9 for HAM since later bursts show rates lower than 
observed. 

The oxygen abundances for these same models are 
shown in Fig. [2] In this plot we can see that the two ef- 
ficiencies chosen, 33 per cent and 10 per cent, give the up- 
per and lower limits respectively for the observed oxygen 
abundance range in this type of galaxies. The models with 
the same initial efficiency and different attenuation type, are 
very similar, and if At is shorter, the low oxygen abundance 
limit is reached very quickly. 

Figs. [3] and [4] show diagnostic diagrams involving the 
ratios of intense emission lines. The number of the model 
(see Table 1) is given in each panel. The higher the effi- 
ciency, the more ionizing photons produced and the higher 
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r(Gyr) 

Figure 2. Evolution of oxygen abundance for same models and 
with the same line code as in Fig. [T] The dashed (green) lines 
define the observational data range 



ionization parameter, leading to a higher excitation of the 
gas. The low efficiency models cover the region of the 
diagram occupied by the young and less metallic galax- 
ies (high [OIII]AA5007,4959/H/3 , low [OII]A3727/H/3 and 
low [NII]A6584/Ha), as expected. The high initial effi- 
ciency models reproduce high excitation and high abun- 
dance galaxies, with high [OIII]AA5007,4959/H/3 and high 
[Nil] A6584/Ha. The ionization degree is driven by the ef- 
ficiency of the bursts, not by the attenuation mode or at- 
tenuation factor. Thus, the models with different attenuation 
modes and the same initial efficiency do not show notewor- 
thy differences, since it is the current burst of star formation 
which produces the observed emission lines, and the under- 
lying population does not affect the ionization parameter. 

The attenuation sets the SFR of the successive bursts 
and determines the contribution of the underlying popula- 
tion. A higher attenuation implies a larger contribution of 
previous bursts to the total SED. Therefore, it is important 
to see the attenuation factor effect on the the evolution of 
the broad-band continuum colors. 

In Fig. [5] we show the evolution of the modeled contin- 
uum colors for each burst of star formation. The inclusion of 
nebular emission continuum in the computation of the SSPs 
reddens the colors of very young populations significantly, 
mainly at low metallicity as explained in MGVB09. 

The HAM cases, compared with SAM, have a higher 
contribution from the non ionizing underlying continuum, 
which makes the colors redder for subsequent bursts and 
the young and blue population features disappear. Mean- 
while, SAM maintains blue colors, characteristic of the cur- 
rent burst of star formation. 

In Fig. [6] the evolution of the colors is shown includ- 
ing the emission lines contribution to the wide band fil- 
ters. We have done so taking into account the strongest 
emission lines that contribute to the color in each broad- 
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Figure 3. Diagnostic diagram for the computed models involv- 
ing [OIII]AA5007,4959/H/3 vs. [OIIJA3727/H/3 of low (top panels) 
and high (bottom panels) initial efficiency models correspond- 
ing to the SAM at the left and to the HAM at the right and 
center panels. The different colored lines represent each burst, 
from the first one occurred at t=0 Gyr (black line) to the last 
one at t=13Gyr in models 1 to 4, at t=1.3Gyr in model 5 and 
t=0.6 Gyr in model 6 (orange line). Observat ional data are from 
iHovos fc Dia3 (2006) —open red squares- and llzotov et af I hood ) 
-grey dots-. The error of the observational data are less than 1% 
of the line intensity, and therefore . they arc not included in the 
figure. 
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log([NIl]/Ha) 
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Figure 4. Diagnostic diagram for the computed models involv- 
ing [OIII]AA5007,4959/H,3 vs. [NII]A6584/Ha. The different col- 
ored lines and dots has the same meaning as Fig [3] The error of 
the observational data are less than 1% of the line intensity, and 
therefore , they are not included in the figure. 



band spectral interval at redshift zero. These are mainly 
[OII]AA3727 in U, H/3 in B, [OIII]AA5007,4959 in V, Ha in 
R and [SIII] AA9069,9532 in /. In that case the more intense 
the last burst, the larger the change in the colors, so the 
dotted lines are more similar to colors calculated only with 
the stellar populations than those shown by solid lines for 
higher efficiencies. On other hand, not all bands are equally 
modified at the same time. 

A fine tuning of the attenuation factor can be obtained 
from the evolution of the broad-band continuum colors ver- 



© 2011 RAS, MNRAS 000, []>?? 



8 Martin-Manjon et al. 




log T 

Figure 5. Evolution of the continuum colors, U-B, B-V, V-R and R-I for each burst along lOMyr. The SAM are plotted on the left, 
the HAM with At = 1.3 Gyr on the intermediate panel and HAM with At = 0.1, and At = 0.05 Gyr, respectively, on the right of the 
diagram. All panels include the low and the high efficiency models as dashed and solid lines, respectively, excepting for the models on the 
right panels, wich are showing both high efficiency models with different inter-burst times, 0.1 and 0.05 Gyr as solid and dashed lines, 
respectively. Each burst is represented with a different colour: from the first burst (black line) to the last one (orange line). 



sus the equivalent width of H/3, EW(H/3). These are im- 
portant observations which give also an effective method to 
uncover the pre s ence of old underlying stellar populations 
|Terlevich et all . 120041 ). 

The comparison of models and observations shows that 
in most star-forming dwarf galaxies, SSPs are unable to 
reproduce the reddest colors shown by the data with the 
lowest EW(H/3) values. In fact the observed H/3 equiva- 
lent width values and colors require a large contribution 



from previous stellar gener ations |Terlevich et all |2004| ; 
iMartfn-Manion et all [200881 ). The evolution of EW(H/3) vs. 
a pseudo-color of the continuum, the intensities of the ad- 
jacent continua of [OII]A3727 and [OIIIJA5007 lines (similar 
to U-V) of our models compared with the data is shown in 
Figure El The SAM and HAM cases for At = 1.3 Gyr are 
plotted in the two first panels. It can be seen that to repro- 
duce the trend of observational data, shifted to red colors at 
low values of EW(H/3) with respect to the SSP predictions, 
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Figure 6. Evolution of the continuum colors, U-B, B-V, V-R and R-I for each burst along lOMyr by including the contribution of the 
emission lines. Each model and each burst are represented with a different colour and line type similar to Figure [5] 



the contribution to the total continuum of the non-ionizing 
red population must be much higher than the contribution of 
the current burst which creates the emission line spectrum. 
This trend can not be reproduced just with SSPs neither by 
increasing their metallicity or age separately nor simultane- 
ously (MMDT). A strong attenuation in the SFR is needed, 
as can be seen in the center panel, to reproduce the whole 
range in EW(H/3) and continuum colors simultaneously. 

The time between bursts (At) is a parameter which 
may also have an effect on the models similar to the attenu- 
ation. The reduction of the time between bursts offsets the 



effect of increasing the attenuation: The underlying popula- 
tion is less evolved and produces less reddening. In Figure [71 
right panel, HAM with At = 0.1 Gyr and At = 0.05 Gyr are 
plotted. The EW(H/3) decreases rapidly while resulting col- 
ors are not shifted to the red sufficiently to cover the range 
shown by the observational sample. The colors of the mod- 
els with shorter inter-burst time are more similar to those 
with soft attenuation (strong bursts), and require an extra 
reddening to reproduce the needed effects of the underlying 
non ionizing populations. However, the EW(H/3) decreases 
more from burst to burst than in the case of a soft attenu- 
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Figure 7. EW(H /3) vs. Iog(l373o/l5oio) compared with observational data from IrTovos fc Di'aj ll2006f ). iTerlevich et al, (1991) and 
ISalzer et al,l (|l995l l . for different attenuation parameters: SAM (left panel), HAM (central panel), and including both, high and low, 
initial efficiencies, shown as solid and dashed lines respectively . The third panel, on the right, corresponds to high efficiency HAM with 
different inter-burst time, At = 0.1 Gyr (solid blue lines) and At = 0.05 Gyr (dashed magenta lines). 



ation, where it maintains a high value at the beginning of 
every bursts. In order to reproduce the observed trend, the 
inter-burst time must be longer than 100 Myr. 

In brief, to reproduce the observable characteristics of 
star-forming galaxies we have to adjust the three input pa- 
rameters: 

• The initial efficiency must be between 0.10 and 0.60 in 
order to produce a first burst which provides oxygen abun- 
dances within the observed values. In fact, efficiencies lower 
than 10% seems more probable to reproduce most of obser- 
vations. 

• Most (perhaps all) Hn galaxies require the contribution 
of previous stellar generations to explain the observed trends 
shown by their continuum and emission line properties. The 
chemical and spectro-photometrical parameters (equivalent 
widths and colors) obtained by our models reproduce the 
observed relations if the contribution of the underlying pop- 
ulation from previous bursts to the total continuum is higher 
than the contribution of the current burst of star formation 
which dominates the observed emission line spectrum. This 
implies a history of star formation higher in the past than 
at present, and, even in that case, attenuated burst along 



the time. The fine tuning of observations is obtained by the 
adjusting of the attenuation factor, once the initial efficiency 
is fixed, (see Martm-Manjon et al. 2012 for details) 

• The inter-burst time must be shorter than 1.3 Gyr and 
longer than 100 Myr. With shorter periods than 100 Myr the 
underlying continuum is too blue and its contribution to the 
color do not reproduce the trend shown by Hn galaxies. 



4 DISCUSSION 

4.1 The different star formation scenarios. 

Three star formation scenario are usually discussed in rela- 
tion with the evolution of dwarf galaxies, (a) Burst: short 
star-formation episodes with large quiescent periods, (b) 
Gasp: long moderate star-formation episodes with short qui- 
escent periods and (c) Continuous or almost continuous star 
formation with few over-imposed sporadic bursts. The three 
cases can be simulated by our models simply with a change 
of parameters: 

• The star-bursting scenario (a) is this one shown in the 
previous sections. The effects of an instantaneous burst can 
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be seen in the emission lines during 10 Myr after the burst 
takes place. Long quiescent periods or, in our case, null star 
formation periods, have been considered to last more than 
1 Gyr, which would be the minimum age of the underlying 
population belonging to the previous burst. Bursts occurred 
before the immediately previous one have very poor contri- 
bution to the current total continuum luminosity, and the 
use of this long time for the inter-burst period produces a 
similar result to assuming a two-burst model for some pa- 
rameters, as EW(H/3). However, differences appear between 
both scenarios since the ISM is enriched by every single 
burst occurring during the galaxy life time, stellar popula- 
tions older than 1 Gyr should be present, and in consequence 
colors would be shifted to the red. 

• The gasping star formation scenario (b) can be sim- 
ulated by substantially increasing the attenuation of the 
burst and reducing the inter-burst time. With the change 
in these two parameters we make a smoother and more 
moderate the star formation. Since the subsequent bursts 
are closer in time, the quiescent periods are shorter. Mod- 
els with extreme attenuation and inter-burst times shorter 
than 100 Myr can be considered a good approximation to 
gasping models. This particular choice of parameters has 
important consequences when compared with the observa- 
tions. Increasing the attenuation a major contribution to 
the continuum by the previous bursts is obtained. How- 
ever, if, in addition, the inter-burst time is reduced, the un- 
derlying population will become younger and bluer. This 
approach would predict the presence of intermediate age 
stellar populations, younger than 1 Gyr, even after several 
star bursts. With larger inter-burst times only the imme- 
diately previous stellar generation contributes substantially 
to the present continuum, but even older stellar generations 
would be noticeable if we further reduce the time between 
bursts. Although complicated, it is also possible to simulate 
more extended star-burst phases by concatenation of sev- 
eral star-bursts. However, population synthesis modeling of 
BCD sp ectra does not favor exten ded periods of star for- 
mation (|Mas-Hesse fc KuntlJ.ll999h . Then, if the starburst 
phase is short, there should be man y of these episodes dur- 
ing th e galaxy lifetime. According to lSanchez Almeida et ahl 
l|2008t l. there should be one BCD phase each 0.3 Gyr, which 
agrees with our models. 

• Under a continuous star formation scenario (c), the SFR 
must be very low and extended. If we reduce the inter-burst 
time to a minimum value and, simultaneously, we reduce the 
intensity of each burst, we obtain a very low star formation 
rate, not comparable to a starburst, but similar to a con- 
tinuous star formation history. In fact, many of the observ- 
able features of star-forming galaxies can be modeled with a 
young stellar population 3-5 Myr old a nd most massive stars 
being essentially coeval. However, as iMas-Hesse fc Kunthl 
l|l999l ) demonstrated, it is also possible to obtain similar 
results with a continuous star formation, lasting at least 
20 Myr since the beginning of the star formation. A low 
continuous star formation rate cannot be neglected, espe- 
cially in low metallicity galaxies. Otherwise, if the SFR 
is very low, (log (SF R) < -3 approximately, according to 
iMartm-Manionl . l2009l .models) . there are not enough ioniz- 
ing photons to produce emission lines. In this case we are not 
reproducing the so-called Hn or BCD galaxies, but another 
type of galaxy or another stage of their evolution. 



4.2 Connecting star-bursting models with the 

different evolutionary phases of a star-forming 
dwarf galaxy. 

Actually, we may consider that, instead of a different sce- 
nario, the low and continuous star formation could be the 
phase occurring between bursts. In this way we could change 
the perspective of the problem and study the characteristics 
of the galaxies for the two phases: the burst and the fol- 
lowing ~ 10 Myr, and the inter-burst, low activity periods, 
that is, after the first 10 Myr from the last star formation 
episode. In fact this is related with the two different time 
scales involved in the evolution, and related with different 
kind of data: On the one hand, a short time-scale in which 
the emission lines and othe effects of the ionizing stellar 
populations are observed. On the other hand there is long 
time-scale which defines the red color evolution mainly due 
to the the stellar populations older than 10 Myr. 

Since a Hn galaxy can be considered just a star- 
forming phase or a stage in the evolution of a gas rich 
dwarf galaxy, we can consider the inter-burst time in our 
models as a quiescent s tar-fo rming dwarf (QBCD) phase 
l|Sanchez Almeida et all l200Sf ) . Then the Hn galaxy may 
be in the phase 10 7 yr after the star burst, in which its ef- 
fects are still visible (emission lines or changes in metallicity 
due to massive stars ejections) and the QBCD phase would 
correspond to the inter-burst periods. The contribution of 
the young and old stellar populations in each galaxy may 
be obtained from the comparison of data coming from the 
two different time-scales, in particular, from observations as 
emission lines, proceeding from the ionizing stellar popula- 
tions, with colors or abundances given by the older than 
10 Myr stars. 

As we have already se e n in ou r previous works, MMDT 
and iMartfn-Manion et.a 3 (|2008bl ), and we also show in 
Fig. |51 when both types of stellar populations are contribut- 
ing to the light, as occurs in BCD galaxies, the colors are 
contaminated by the emission lines, showing in most cases 
different trends than those standard sequences defined by 
the stellar populations, either young or old. Taking into ac- 
count this contamination, the colors observed are those cor- 
responding to the ionizing continuum of a given galaxy, al- 
ways during the BCD phase, that is, the colors of the young 
population which dominates the emitted light. In order to 
study the colors in the inter-burst periods and to know the 
properties of the underlying population, we should further 
eliminate the continuum arising from the starburst and iso- 
late the properties of the underlying host galaxy. 

Figures[UHand[lO]show the B " V vs - V" 1 . V" 1 vs - R " : and 
g-r vs. r-z color-color diagrams for our 6 models. Black cir- 
cles show the stellar continuum colors in all panels, while the 
colors computed including the contribution by the strongest 
emission lines are shown as green and red for models 1 and 
2 respectively, brown and orange for models 3 and 4, and 
blue and magenta for models 5 and 6. In order to repre- 
sent these older phases we have also represented, as colored 
open squares, the colors corresponding to SSPs, as given in 
MGVB09, for ages older than 10 Myr that is, without tak- 
ing into account the youngest stellar populations phases, and 
metallicities from Z=0.0001 to 0.02. We have not computed 
exactly the colors of our inter-burst periods, but they must 
fall in this same region of the diagrams. 
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Figure 8. B-V vs. V-I color-color diagram for our models including the stellar continuum colors of the starburst phase (black dots), and 
these continuum colors contaminated by the emission lines. In the top panel we have green and red dots for low and high efficiency SAM 
models. In the middle panel, models 3 and 4 as brown and orange dots; And in bottom panel models 5 and 6 with blue and magenta 
dots. Colors of the inter-bur s t per iods (age > 10 Myr) for 6 metallicities are the open (colored) squares. Observational data arc taken 
from lSanchez Almeida et al] l|2008h -grey small dots- and from lCairos et a l. (2002, 2001a) -cyan full dots and green open and fill triangles 
(see text for an explanation). 



In the B-V vs. V-I diagram (Fig. [S]) the inclusion 
of the emission lines contribution to the continuum colors 
shifts the position of the model points almost perpendicu- 
larly to the originally computed ones. The location of the 
points is mainly determined by the contribution of strong 



[OIII] AA5007,4959 emission lines to the V band, and, logi- 
cally, they are closer to the normal stellar populations locus 
when efficiency is low. This means that the excursion from 
locus defined by the main sequence for stellar populations is 
stronger when the star formation efficiency is higher, lead- 
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ins the poin t s far. W e also show in this figure the data by 
ICairos et al.l (|2001al lbT) as cyan dots with error bars. This 
set of observations corresponds to readily observed colors, 
including both continuum and line emission, and no redden- 
ing correction has been applied. It can be seen that some 
of the data are impossible to be reproduced by the models 
which do not take into account the contribution by emission 
lines, even if some am ount of reddening is invoked. On the 
other hand the data bv lCairos et al.l (|2002l ) of resolved loca- 
tions in Mrk 370, shown in the figure as open green triangles 
are very well reproduced by our continuum colors. The data 
of Mrk 370 are based on UBVRI broadband and Ha nar- 
row band observations. In this case, the authors subtracted 
the contribution of the underlying continuum from the old 
stellar population, removing the contribution from emission 
lines and correcting for extinction, measuring in this way 
true colors of the young star- forming knots. This way while 
open green triangles correspond to the observed Mrk 370 
colors, including the emission lines contribution and the un- 
derlying stellar population, the full green triangles are the 
stellar population colors after correction of extinction, emis- 
sion lines and old stellar populations. In this case we see 
that some points fall out of the expected region of models 
which we interpret as probably due to a over-correction with 
SSPs models. The underlying component has colors redder 
than the BCD ones, and they are also well reproduced by 
our models. 

In Fig. [9] we pl ot V-I vs. V -R, w ith the same 
data and code from ICairos et"aH (|2001al lbl. 120021 ') and 
ISanchez Almeida et al.l ( 20081 ') than in the previous figure. 
Moreover, we have also included th e data of BCDs and its 
hosts from lTelles fc Terlevichl (|l997l V This sample consists of 
15 BCDs for which they observed the total colors. The data 
corrected of extinction are the green stars. When they cor- 
rect of the emission lines effect, the obtain the colors shown 
by magenta stars. These points corresponding to the ac- 
tual young stellar populations are closer to the region where 
our models are represented by the same coloured dots as 
in Fig. [S] They subtracted the contribution of the underly- 
ing galaxy by assuming the mean surface brightness of the 
extensions, or hosts, to be constant and representative of 
the underlying galaxy within the starburst regions. Sum- 
marizing, they give separately the colors of this underlying 
stellar continuum (blue stars) within the starburst, the con- 
tributions including the emission lines (green stars) and the 
colors of the BCD stellar continuum (magenta stars), cor- 
rected from emission lines, which are close to black dots, 
corresponding to our continuum BCD phase. The colors of 
the isolated host galaxy (blue stars) show redder colors than 
the BCD component, in agreement with QBCDs, Mrk 370 
host and model colors with ages older than 10 Myr. Again 
we note some differences with our models that we assign to 
an over-correction of the emission lines. 

Finally we represent in Fig llOl the g-r vs. r-z diagram 
with the same symbols that before. The selected QBCDs 
(|Sanchez Almeida et al. . 20081') are assumed to be like BCD 



host galaxies (|Amorin et al. 020071 ; I Amorm et al.l . [2009l ). The 

colors of the SSPs with ages older than 10 Myr (colored 
squares) can reproduce these observations, indicating that 
they correspond to the colors of a more evolved population 
than the one corresponding to the current burst in BCDs. 
However, some of the observed dots lie in the BCD zone, 



indicating that this data sample also contains a certain pro- 
portion of BCDs, following a continuous evolutionary se- 
quence. In this color-color diagram there is a larger differ- 
ence than in the previous ones between models and data for 
the bluest points, showing these data higher g-r values than 
predicted for similar low r-z colors. Probably this is due to 
the transformation used to calculate SDSS colors with our 
Johnson colors, which would be not valid for this type of ob- 
jects with emission lines contaminating the wide band filter 
magnitudes. 

We see in the three figures that models start to separate 
off the main sequence of colors at different points following 
the type of attenuation. In the top panel, there are points 
with large contribution of lines starting even at the lowest 
and left corner, where the youngest stellar populations lie. 
The same points in medium panel begin to separate at red- 
der colors, this implies that the underlying continuum corre- 
sponds to an older age. These models cover all observational 
range. However in the bottom panel the points with the con- 
tribution of the lines fall in a region much bluer where there 
are no observations. The underlying stellar populations are 
much younger than the populations of the data. 

In any case it seems that BCDs and QBCDs colors 
are well reproduced by the same basic models, during the 
first 10 Myr after the beginning of star formation and after 
10 Myr up to more than 1 Gyr respectively, and therefore we 
can say that both types of galaxy colors overlap and could 
be considered as different phases in the evolution of the same 
object. 



5 SUMMARY AND CONCLUSIONS 

Historically, the most common scenario assumed for the star 
formation history of the star-forming galaxies is that stars 
form in a bursting mode. This scenario consists in successive 
instantaneous star formation bursts interposed between long 
periods of null or very low star formation activity. If these 
periods are shorter and the bursts are more extended and 
moderate, we are under a gasping star formation scenario. 
On the contrary, if the star formation is low and moderate 
during the whole life of the galaxy, and only some sporadic 
bursts take place, a continuous star formation scenario is 
appealed. 

We have analyzed the possibilities of these scenarios by 
means of theoretical models based on the successive bursts 
star formation hypothesis. Each galaxy is modeled assuming 
an initial amount of unprocessed gas of 10 s Mq . The evolu- 
tion is computed along a total duration of 13.2 Gyr during 
which successive star-bursts take place. Since any of the pos- 
sible scenarios of the star formation histories described show 
evidences of being able to explain the observed data, we have 
tried to check the hypotheses suggested. 

The grid of theoretical models are computed by the 
combination of three tools: a chemical evolution code, an 
evolutionary synthesis code and a photo-ionization code, all 
of them previously calibrated. They have been used in a self- 
consistent way, i. e. taking the same assumptions about stel- 
lar evolution and nucleosynthesis, and the resulting metal- 
licity in every time step. These models have three free input 
parameters which can be changed to obtain different model 
results: (a) The initial star formation efficiency which deter- 
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Figure 9. V-I vs. V-R color-color diagram for our models including the stellar continuum colors (black dots), and the continuum colors 
contaminated by the emission lines. In the top panel we have green and red dots for low and high efficiency SAM models. In the middle 
panel, models 3 and 4 as brown and orange dots; And in bottom panel models 5 and 6 with blue and magenta d ots. Colors of the inter- 
burst periods (age > 10 Myr) for 6 metallicities a re the open ( colored) squares. Observational data are taken from lSanchez Almeida et al 



200S ) -grey small dots- and from ICairos et al ] l|200l 12001 J ) -cyan full dots and green open and fill triangles . and lTel!eT^r*^rlevich1 



19971 ) -green and magenta stars— ( see text for explanations). 



mines the initial star formation rate (SFR) and the initial 
metallicity of the gas, (b) the attenuation of the successive 
bursts, which determines the evolution of the gas, keeping 
metallicity between the observational data limits, and (c) 



the inter-burst time, which sets the age of the non ionizing 
underlying population. 

We have shown that (a) also leads the behavior of the 
ionizing gas since emission lines are determined by the num- 
ber of massive stars and the metallicity of the emitting gas. 
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With this parameter we can control the number of old pris- 
tine stars that we are going to find in the next star formation 
burst. 

The strength of the bursts is controlled by (b). The more 
attenuated, the less the star formation rate, which favors the 
previous stellar generations to have more weight on the total 
current spectrum, hence modifying the equivalent width and 
the colors of the continuum. 

Finally, the contribution of the underlying population 



can be determined by (c). These periods are the quiescent 
phases of the BCD. The shorter the inter-burst time, the 
younger the underlying stellar populations. These stars will 
be part of the host BCD galaxy, that is, all those generations 
of stars which have not ionizing stars, showing ages of more 
than 10 7 yr. 

The models can reproduce simultaneously the data rel- 
ative to the current ionizing population and the data which 
give us information about the evolutionary history of the 
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star formation. In order to reproduce the observational data, 
the parameters are constrained: 

a) the initial star formation efficiency must be lower than 
50% and more probably is around 20% . 

b) the subsequent star-forming burst must be attenuated, 
in order to have the adequate proportion of underlying to 
young/ionizing stellar populations; and 

c) the inter-burst time must be longer than 100 Myr but 
shorter than 1.3 Gyr. 

With our models we are able to compute the stellar 
continuum colors as well as the broad band colors corrected 
from the contribution of the more intense emission lines. 
This effect shifts the colors almost perpendicularly to the 
stellar continuum ones, specially during the first bursts of 
star formation, when the emission lines are more intense. 
The contamination of the continuum by the emission lines 
can make the observed galaxies appear younger than they 
really are, hiding evidence of underlying old populations. 

On the other hand, the different characteristics shown 
by the star formation histories of our models, imply that we 
may be observing the same objects in different evolutionary 
stages, that is, just during the burst, in the post-burst phase 
or in the inter-burst periods. If we consider that the differ- 
ent star formation scenarios correspond to different phases 
in the life of a single galaxy, our models are able to reproduce 
the characteristics of each one of these phases. Considering 
the BCD phase as the first 10 Myr after a star formation 
burst, and the quiescent phase as the inter-burst periods, 
we can reproduce the galaxies both in a high activity and 
in a post-starburst (or pre-starburst) phase. The main dif- 
ference between these phases are the existence of emission 
lines produced by the ionized gas. Their intensity, and their 
contribution to the total continuum color, can be modified 
by the initial efficiency. With the attenuation factor we can 
produce a more or less reddened stellar continuum due to 
the contribution of the host galaxy, the age of these old stars 
being determined by the inter-burst time. 
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